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An analytical expression of probability density function 
(PDF) of velocity fluctuation is derived with the help of the 
statistics based on generalized entropy (the Tsallis entropy or 
the Renyi entropy). It is revealed that the derived PDF ex- 
plains the detailed structure of experimentally observed PDF 
as well as the scaling exponents of velocity structure func- 
tion. Every parameters appeared in the analysis, including 
the index proper to the Tsallis entropy or the Renyi entropy, 
are determined, self-consistently, by making use of observed 
value of intermittency exponent. The experiments conducted 
by Lewis and Swinney (1999) are analyzed. 
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An investigation of turbulence based on the general- 
ized entropy, Tsallis' p]|| or Renyi's ||, was started by 
the present authors pF with an investigation of the p- 
model HQ in terms of a generalized statistics constructed 
on the entropy. We intend to analyze the velocity fluctu- 
ation, 5u(r) — \u(x + r) —u(x)\, between two points with 
a distance r apart in a turbulence. Here, u represents 
one of the components of the fluid velocity field u, say 
x-component. By making use of the velocity fluctuation 
5uq of eddies with the largest size £q, the Reynolds num- 
ber is estimated, in the absence of intermittency M, as 
Re = 8uq£q/v = (A)/«) 4 ^ 3 where v and r\ — (v z / 'e) 1 ™ are, 
respectively, the kinematic viscosity and the Kolmogorov 
scale. The quantity e represents the energy input rate at 
the largest eddies. 

In the case of high Reynolds number Re > 1, there 
exist a lot of steps, n = 1, 2, • • •, at each of which eddies 
break up into two parts producing a energy cascade. The 
size of eddies in the nth step of the cascade is assumed 
to be given by £„ = S n £o with S n — 2~™. Then, our 
main interest in the following reduces to the fluctuation of 
velocity difference 8u n = Su(£ n ) corresponding to the size 
of nth eddies in the cascade. Note that the dependence of 
the number of steps «k on r/rj, within the analysis where 
intermittency is not taken into account M, is given by 



n = - log 2 r/rj + (3/4) log 2 Re. 



(1) 



In this paper, we examine the experimental results ob- 
tained by Lewis and Swinney Jl6[ for turbulent Couette- 
Taylor flow produced in a concentric cylinder system. 



Our analysis is based on the fact that, for high Reynolds 
number Re 3> 1, the Navier-Stokes equation for in- 
compressible fluid is invariant under the scale transfor- 
mation M: r — > A Q / 3 it, u — > Xr, t — > A 1_a / 3 t and 
(p/p) ~ * A 2 "/ 3 (p/p)- Here, the exponent a is an arbi- 
trary real quantity which specifies the degrees of singular- 
ity in the velocity gradient \du(x) / dx\ = lim^ n _,o \u(x + 
£ n ) — u(x)\/£ n = lim£ n ^o Su n /£ n . This can be seen with 
the relation 5u n /5uo — (£ n / %) Q ^ 3 ', which leads to the 
singularity in the velocity gradient || for a < 3, since 

8u n /£n^ — const.. 

The distribution P^(a)da oc P^(a) n da of singular- 
ities at the nth step in the cascade with 



P«(a)oc [1 - (a - a ) 2 /(Aa) 



21 1/(1-9) 



(2) 




(Aa) 2 = 2X/[(1 - g)ln2] was derived 
taking an extreme of the Tsallis entropy 
Sj[PW(a)} = (1 -q)' 1 (J daPW (a)" - l) which is 
non-extensive, or the Renyi entropy || S^[P^ (a)] — 

(1 — q) 1 In J daP( 1 \a) q which has the extensive charac- 
ter, with appropriate constraints, i.e., the normalization 
of distribution function, JdaP^(a) = const., and the 
q- variance being kept constant as a known quantity, a 2 q = 

((a-a ) 2 ) q = {JdaPW(a)i(a-a ) 2 )/ fdaPV>(a)«. In 
deriving PW(a), we assumed that each step in the cas- 
cade is statistically independent. Note that the values 
of a are restricted within the range [a m i n , a max ], where 
ctmax — a = aa - a miu = Aa. 



fit) 



we 



With the help of the relation P^(a) oc Sn 
can extract the multi-fractal spectrum /(a) in the form 



/(a) = 1 + (1 



q) 1 iog 2 



1 - (a - a ) 2 I (Aa) 2 , (3) 



which reveals how dense each singularity, labeled by a, 
fills physical space. 

By making use of an observed value of the intermit- 
tency exponent fi as an input, the quantities ao, X and 
the index q can be determined, self-consistently, with the 
help of the three independent equations, i.e., the energy 
conservation: (e n ) = e, the definition of the intermit- 
tency exponent u: (e n ) = e 2 8~ tl , and the scaling rela- 
tion [|-|^): 1/(1 - q) = 1/a- - l/a+ with a± satis- 
fying f(a±) = 0, where the average (■ • ■) is taken by 



1 



p( n \a). The scaling relation is a generalization of the 
one derived first by 0,^5) to the case where the multi- 
fractal spectrum has negative values. The /^-dependences 
of the self-consistent solutions, ao, X and q, are given 
in Fig. |l] and Fig. || for the region where the value 
of \i is usually observed, i.e., 0.16 < /x < 0.31. We 
see that they are given in this region by the equations 
a = 0.998 + 0.587/x, X = -5.73 x 10~ 3 + 1.21// and 
q = -0.607 + 6.01^ - 7.72/r 2 . 
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FIG. 1. The //-dependence of ao and X. Squares and cir- 
cles are the points where the self-consistent equations are 
solved. 




FIG. 2. The //-dependence of q. Circles are the points 
where the self-consistent equations are solved. 

Adopting the value /i = 0.28 observed in the ex- 
periment |ll| with the Taylor-scale Reynolds number 
Rx = 262 (Re = 540 000) or R x = 80 (Re = 69 000), 
and solving the above three equations self-consistently, 
we have q = 0.471, a = 1.162, X = 0.334. Then, we 
obtain a+ - a Q = a Q - a_ = 0.748, Aa = 1.350 

It may be reasonable to assume that the probability 
II^ n > (xn)dx n to find the scaled velocity fluctuation \x n \ = 
Sun/Suo in the range x n ~ x n + dx n can be divided into 
two parts: 

n^\ Xn )dx n = n^\x n )dx n + n^\\x n \)dx n (4) 



where the singular part Z7g (|x„|), stemmed from mul- 
tifractal distribution of singularities in velocity gradi- 
ent that may be related to the fluctuations caused by 
turbulent viscosity, is determined by i7g (|x„|)dx n = 
p( n >(a)da with the transformation of the variables be- 
tween \x n \ and a: \x n \ = 5"^ 3 , leading to Q 

i7 s (n) (|x„|)dx n = 3G^(\x n \)/ (zW\]n5 n \) dx n (5) 

with 

G^{x) = x- 1 [1 - (a(x) - a ) 2 /(Aa) 2 ] n/il - q) , (6) 

a(x) = 31nx/lnc)„. Note that the singular part of the 
probability density function has a dependence on In 5u n , 
and that the values of I restricted within the range 

[<5" max ^ 3 , (5" min ^ 3 ]. The first term, II^ l \x n ) assumed to 
come from thermal or dissipative fluctuation caused by 
the kinematic viscosity, will be considered later in this 
paper after estimating the moments of the velocity fluc- 
tuation. 




FIG. 3. The scaling exponent f m of the velocity structure 
function. The experimental results obtained by Lewis and 
Swinney are shown by circles (Re = 540 000) and by triangles 
(Re = 69 000). The present theoretical result ([]) is drawn by 
solid line for the intermittency exponent /i = 0.28 taken from 
the experimental data. Dotted line represents K41, whereas 
dashed line She-Leveque. The prediction of the log-normal 
model is given by short-dashed line with the best fit but incon- 
sistent value jj, = 0.27 of the intermittency exponent following 
Lewis and Swinney. 

The mth moments of the velocity fluctuation, defined 

by ((\x n \ m )) = JZodXnK^n^ixn), are given by 

((Kn)=2 7 W + (l-2 7 (" ) )a m e (7) 

where a M = {2/[C\ /2 {l + C^ 2 )]} 1 / 2 with C q = I + 
2q 2 (l - g)Xfn2, and 

2 7 W = / dx n \x n \ m Il£\x n ). (8) 
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Here, we used the normalization ((1)) = 1. The quantity 
C m = a m/3 - 2Xm 2 / [9 (l + C%fy 



with 77 s (n) = 3(1 - 2^ n> )/(2i ny /2nX\]nS n \). 

In the following, we will derive the symmetric part of 
the PDF of velocity difference. The origin of the asym- 
metry in the PDF (i.e., skewness) may be attributed to 
dissipative evolution of eddies or to experimental setup 
and situations. The consideration of the latter will be 
reported elsewhere. 



1 - log 2 ll + C 



,1/2 
y m/3 



/(I -9), 



(9) 



is the so-called scaling exponents of velocity structure 
function, whose expression was derived first by the 
present authors P-jl^l- In Fig. [|, we compare the present 
result (|^) with the experimentally measured scaling ex- 
ponents Jl|] at Re = 540 000 (circles) and Re = 69 000 
(triangles). We use the observed value of the intermit- 
tency exponent, i.e., fi = 2 — ( & = 0.28, for theoreti- 
cal analysis. There are also represented, as references, 
the predictions of K41 (dotted line) 0, of She-Leveque 
(dashed line) |l7[] and of the log-normal (short-dashed 
line) |lg| -|2C|] . Note that our formula (^|) can also explain 
the scaling exponent for higher moments M Q . 




FIG. 4. The n-dependence of the PDF ij (n) (£ n ) given by 
the present analysis with q = 0.471 (/1 = 0.280) for integer 
values n from 5 to 20, from left to right at the axis £. 



Let us introduce new variable 



scaled by the variance of velocity fluctuation where = 
[2j^ n) S-^ + (1 - 2^ n) )a 2 }- 1/2 , and the probability den- 
sity function (PDF) -ff < -™- ) (l£« I) in terms of this variable 
through the relation 



(11) 



Making use of the relation between and a, the PDF 
responsible for the distribution of singularities TlS (|£n|) 
can be rewritten in terms of a in the form 




FIG. 5. Experimentally measured PDF of the velocity fluc- 
tuations by Lewis and Swinney for _Ra = 262 (Re = 540 000) 
are compared with the present theoretical results 77 ( '™'(£n). 
Open triangles are the experimental data points on the left 
hand side of the PDF taken from Lewis and Swinney, whereas 
open squares are those on the right hand side. Closed circles 
are the symmetrized points obtained by taking averages of 
the left and the right hand experimental data. For the ex- 
perimental data, the distances rjr\ = t-n/rj are, from top to 
bottom: 11.6, 23.1, 46.2, 92.5 208, 399, 830 and 1440. Solid 
lines represent the curves given by the present theory with 
q = 0.471. For the theoretical curves, the number of steps in 
the cascade n are, from top to bottom: 14, 13, 11, 10, 9.0, 
8.0, 7.5, 7.0. For better visibility, each PDF is shifted by -1 
unit along the vertical axis. 



In order to determine the normal part ilj^ (£ n ) of the 
PDF, we divide 77 (n) (Cn) into two parts |HJ: ffW (£„) = 
n [ <KU) for \U < C, and = for 



n /3-C 2 /2 



The point £* is defined by 



(n) 



where a* is 



77 S M(|^|) = fl£OjC a /2-a/3+W(aO 



(12) 



a < |& I < inSn" 

C=6X* /3 - C2/2 giving 77 s (n) (|Cl) 
the solution of C2/2 - a/3 + 1 - f(a) = 0. The PDF func- 
tion 77<* (£n) is constituted both by thermal fluctuation 
and by the multifractal distribution of singularities, and 
is assumed to be given by Gaussian function On 
the other hand, the contribution to fl^ (£„) is assumed 
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to come only from the multifractal distribution of sin- 



gularities, i.e., 11^ (£ n ) = /7g" ; (|£„|) [|2lJ. Connecting 



(n), 



n. 



(«), 



<<*'(£n) and (£„) at £* having the same value and 
the same derivative there, we have 



[l+3/'(«*)][«»/C) 3 -l]/2 



(13) 



The n-dependence of i7(™)(£ n ) is shown in Fig. |J with 
the self-consistently determined parameters ao, X and q 
pl[ . We notice that there are two points which look like 
independent of the number n of steps in the cascade at 

w 0.8 and at £„ « 2.3. These points are also observed 
in the experimental data Jul at about the same values 
of £„. 

The comparison of experimentally measured PDF's of 
the velocity fluctuation |l6| and those obtained by the 
present analysis is given in Fig. 0. In order to extract 
the symmetrical part from experimental data, we took 
mean average of those on the left hand side (represented 
by open triangles in the figure) and the right hand side 
(by open squares). The symmetrized data are described 
by closed circles. The solid lines are the curves n^ n '(^ n ) 
obtained by the present analysis. 




FIG. 6. Dependence of n on r /rj, extracted from Fig. |^, is 
given by circles, which can be fitted by solid line jl4|). The 
dashed line represents (|l]), put for reference, which is derived 
in the absence of intermittency. 



The dependence of n on r/rj, extracted from Fig. |^, is 
shown in Fig. @ by solid line, which gives us the best fit 



-1.019 x log 2 r/r] + 0.901 x log 2 Re 



(14) 



with Re = 540 000. The dashed line represents (|j]) valid 
in the absence of intermittency. From (|l4]), we see that 
the power w in Re = (io/v) w turns out to be w = 1.13 
in the present intermittent turbulence, which is 4/3 in 
the absence of intermittency, and that the number tik 
corresponding to the Kolmogorov scale is estimated as 
?ik = 17.2 for the experiments under consideration. 



The success of the present theory in the analysis of the 
experiments [|l6| may indicate the robustness of singu- 
larities even for the case of no inertial range. The same 
experiments are investigated in |22j ] by a rather different 
analysis from the present one. Comparison of these two 
approaches is given in pl| . The present theory works 
quite well also for another systematic numerical experi- 
ments conducted by Gotoh [23] ]. It will be reported else- 
where. 

The authors would like to thank Prof. C. Tsallis for 
his fruitful comments with encouragement. The authors 
are grateful to Prof. T. Gotoh for his kindness to show 
his data prior to publication. 
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